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We investigate the critical scaling of the quark propagator of Nf = 2 QCD close to the chiral 
phase transition at finite temperature. We argue that it is mandatory to take into account the back- 
reaction effects of pions and the sigma onto the quark to observe critical behavior beyond mean 
field. On condition of self-consistency of the quark Dyson-Schwinger equation we extract the scaling 
behavior for the quark propagator analytically. Crucial in this respect is the correct pion dispersion 
relation when the critical temperature is approached from below. Our results are consistent with the 
known relations for the quark condensate and the pion decay constant from universality. We verify 
the analytical findings also numerically assuming the critical dispersion relation for the Goldstone 
bosons. 

PACS numbers: 



I. INTRODUCTION 

The phase diagram of strongly interacting quark mat- 
ter has received a lot of attention, both from theoretical 
and experimental perspectives. One of the debated is- 
sues in this respect is the nature of the phase transition 
in Nf = 2 QCD in the chiral limit. It has long been 
suggested that the SU(2) V x SU{2)a — 0(4) symmetry 
in this case is broken spontaneously to SU(2)v — 0(3) 
by a second order chiral phase transition of the 0(4)- 
universality class in 3-dimensions The underlying 
symmetry breaking pattern results from the assumption 
that the anomalously broken Ua(1) symmetry of QCD 
leads to a relatively heavy pseudoscalar flavor singlet me- 
son also for temperatures around T c . 

One of the tools to investigate this issue is lattice gauge 
theory, see e.g. @, The problem is delicate since the 
analysis of a potential second order phase transition in a 
finite volume is elaborate and light quarks are expensive 
in lattice computations. As a result there are indications 
for 0(4) scaling from simulations with Wilson fermions 
[3| and very recently from results with staggered fermion 
actions [5|-|7|. However, since there are also arguments in 
favor of a first order transition || Q , the problem seems 
not yet completely settled. 

The symmetry breaking pattern with an anomalously 
broken Ua{^) symmetry is a widely used assumption for 
the construction of effective models as for example the 
Polyakov-loop quark meson model (PQM), see e.g. [To|- 
[lH . Recent constructions of effective models include and 
investigate the effects of the anomaly fl4l uM ■ T he issue 
of critical scaling has been addressed in fla . fl7| . 

A third approach to the phase diagram are functional 
methods like the functional renormalization group 0- 
l20j or Dyson-Schwinger equations (DSEs) pl| - |23j applied 
directly to QCD. In principle, these methods are capable 
to link information on the microscopic degrees of freedom 



of the theory to the characteristics of the phase transi- 
tions. Recent works include the calculation of the non- 
perturbative Polyakov-loop potential from the propaga- 
tors of Yang-Mills theory []jj, the investigation of the 
chiral and deconfinement transition [20l . u2\ , [23| . quark 
spectral functions (24| and the exploration of color su- 
perconducting phases j2o| . 

Although DSEs are well established they are usually 
not used to study critical phenomena and it may be even 
unclear to what extend they are suitable. In fact, up 
to now critical scaling in DSEs has only been observed 
on the mean field level [26|, [27} • The main goal of this 
paper is to go beyond these studies. We employ the 
Dyson-Schwinger equations to explore the 0(4) scaling 
in the quark propagator at the chiral phase transition 
for Nf = 2 QCD. By including the effects of pion and 
sigma back-reaction onto the quarks along the lines of 
the zero temperature studies in Refs. [HI, [29[ we identify 
a mechanism that generates critical scaling beyond the 
mean field level. Our results therefore provide a basis 
for further studies of the chiral phase transition in this 
framework. They also serve to study the connection be- 
tween the scaling properties in the effective theory and 
the behavior of the fundamental microscopic degrees of 
freedom such as the quark propagator. 

The paper is structured as follows: In section[TT]we dis- 
cuss our approximation scheme for the DSE of the quark 
propagator at finite temperature. We then analyze the 
scaling behavior of the DSEs analytically in section IIIII 
and demonstrate that the mechanism for the onset of 
critical scaling beyond the mean field level in the quark 
propagator is located in the pion and sigma dispersion 
relations. For our numerical investigations we discuss a 
tractable truncation scheme in section IIV Al and verify 
our analytical findings in section llV Bl under the assump- 
tion of a scaling law for the pion velocity. We conclude 
and summarize in section [V] 
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II. QUARK DYSON-SCHWINGER EQUATION 
AND EFFECTS FROM GOLDSTONE BOSONS 

A. Backcoupling mesons onto quarks 

The DSE for the quark propagator is given diagram- 
matically in Fig. [TJ The unknown quantities in this 
equation are the dressed gluon propagator and the one- 
particle irreducible (1PI) quark- gluon vertex. We are 
particularly interested in the behavior of the quarks close 
to the chiral limit and in the vicinity of the chiral phase 
transition. There the relevant degrees of freedom are 
associated with long-range correlations and can be iden- 
tified as the Goldstone modes and the associated radial 
excitation. These contribute to the DSE of the quark- 




FIG. 1: The Schwinger-Dyson equation for the fully dressed 
quark propagator 

gluon vertex as outlined in Ref. J28[. We first briefly 
summarize the relevant findings of Ref. 28[ and subse- 
quently consider their implications at finite temperature. 

The DSE of the 1PI quark-gluon vertex can be an- 
alyzed in terms of contributions from pure Yang-Mills 
theory and unquenching effects in the form of the back- 
coupling of meson degrees of freedom onto the quark- 
gluon interaction [28|, [2^1 ■ To lowest order in a skeleton 
expansion these appear as one-meson exchange between 
quark and anti-quark as displayed in the upper panel 
of Fig. [2] Since we are primarily interested in long-range 
meson contributions, we need to keep this term explicitly 
in our approximation scheme of the quark-gluon vertex. 

At zero temperature and at physical quark masses the 
meson back-reaction term is dominated by the exchange 
of the lightest mesons, i.e. pions. At finite temperature 
and close to the chiral limit the 5'f7(2) y 4-symmetry gets 
restored at the critical temperature T c and the order pa- 



rameter correlator needs to be taken into account. This 
is done by considering the flavor singlet sigma meson as 
well. The pseudoscalar flavor singlet meson, however, is 
assumed to receive a topological mass due to anomalous 
non-conservation of the J7yi(l)-symmetry and is therefore 
suppressed. Consequently, we approximate the quark- 
gluon vertex by a sum of purely gluonic contributions 
not affected by critical modes and the meson exchange 
contribution containing pions and the sigma meson. The 
resulting DSE for the quark propagator is shown in the 
lower panel of Fig. [5J 

The calculation of the quark self-energy with meson 
exchange is in general quite demanding. It contains a 
two-loop diagram containing the full Bcthc-Salpcter am- 
plitudes, Ttt o-, which can in principle be determined from 
the corresponding Bethe-Salpeter equation. In order to 
obtain a tractable truncation that preserves the long 




FIG. 2: The Dyson-Schwinger equation for the quark propa- 
gator using a quark-gluon vertex that includes back-coupling 
effects from pions and the sigma meson. 



range effects close the chiral phase transition we will use 
some low-energy properties of pions in the next section. 



B. The DSEs at finite temperature 

Before wc go into the details, we first give the formal 
expression for the DSE of the quark propagator as shown 
in Fig. [2j At finite temperature it is given by 



S 1 (p,Up) = Z 2 S 1 (p,Up) + T,YU(P, Wp) + S,r(P) + £<r(p, Wp) (1) 

where the self-energies are 

Sym(p,w p ) = Z 1F g 2 C F T^2 J T^yj lv S (q,u q )T™(q,Lj q ,p,u} p ) D^(k,u> k ) (2) 

S w (p,w p ) = 2Z 1F g 2 C F T^ [ ^^^S{q,(J q )Tl(q,u q ,p,(j p )D^(k,u} k ) (3) 

n q 

= Z 1F g 2 C F T^2 J l^S{q,u q )Tl{q,u qi p,uj p ) D^{k,u: k ) (4) 
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with k = p — q = (fc,Wfc), the temperature T, the Casimir 
Cf = (Nc — 1)/(2N C ) and the renormalization factors 
Z\f of the quark gluon vertex and Z2 of the quark prop- 
agator. We used the abbreviations Yv for the triangle 
meson exchange that is part of the quark-gluon vertex. 
This should not be confused with the Bcthc-Salpcter am- 
plitudes r w cr . All color and flavor traces have been car- 
ried out already 1 . In this work we only consider Nf = 2 
and N c = 3. The gluon propagator is denoted by D^ v 
and the purely gluonic part of the dressed quark-gluon 
vertex by T YM . Using a decomposition in Dirac compo- 
nents the dressed quark propagator can be parametrized 
by 



S(p, u> p ) = [-i~f 4 uj p C{p, u> p ) - tfp A(p, u p ) 
= ijiUJp a c (f, io p ) + ijp a a (p, u P ) 



B(p,uj p )} 
-a B (p,u p ) 



(5) 



with the scalar dressing function B(p,uj p ) and the two 
vector dressing functions A(p,ui p ) and C{p,uJ v ). Its bare 
counterpart reads Sq 1 (p,lo p ) = —ij/ + m and we work 
with very small quark masses m — > 0. 

The details of the non-perturbative gluon propagator 
and Yang-Mills quark-gluon vertex are not important for 
the scaling behavior in the vicinity of the chiral phase 
transition. This may be expected from the universality 
argument and will also show up in our analysis presented 
in the next section. In this section we only state the gen- 
eral tensor decomposition of the gluon propagator and 
Yang-Mills quark-gluon vertex that is used in this paper. 
For further details we refer to section HV Al where we con- 
sider the numerical solution. Throughout this paper we 
will work in Landau gauge. In this case the most general 
representation of the gluon propagator at finite temper- 
ature is composed of two dressing functions. It can be 
written as 



D tll ,(k,ui k ) 



Z T (k,UJ k ) T Z L (k,u} k ) t 



'Hvy - uv 

with transverse and longitudinal projectors (i, j = 1 . . . 3) 



k 2 



P£(fc,w*) = P^ikw) - P£,(k,u k ) . (7) 

and dressing functions Zx(k,uj k ) and Zi,(k,u!k)- For the 
Yang-Mills part of the quark-gluon vertex we adopt the 
rainbow ladder approximation 



iYM 



(q,cj g ,p,(jp) ->■ ^ v T YM (k,(J k ) 



(8) 



1 Note that the flavor trace for the self-energy part due to the pion 
backreaction results in a factor of two, since the flavor of the 
external quark is fixed. This has been overlooked in Refs. l28l 
|29|| . where a factor of three has been used leading to a slight 
overestimation of the pion backcoupling effects. 



where r YM (fc, uj k ) is a vertex model depending only on 
the gluon momentum. 

Now we need to specify the meson contributions T^' 7 
in Eqs.© and (HJ). These are given by 

X^ v S{p + f,ui p + u) r )Y va {p,uip,r,ui r )D n ^(f,ui r ) 

(9) 

with the Bcthe-Salpeter amplitudes r WjCT and the meson 
propagators D„ a . In the following we consider the pion 
Bcthe-Salpeter amplitude and propagator. 

In Ref. [3(| Son and Stephanov showed that the real 
part of the pion dispersion relation in the symmetry bro- 
ken phase is given by 



u 2 = u 2 (f + ml) 



(10) 



where u denotes the pion velocity and the pion screen- 
ing mass. At p = the energy of a pion E n = u m n is 
called the pion pole mass. Furthermore as already noted 
in Ref. |3l| there are two distinct pion decay constants 
in the medium. At finite temperature these are given by 
f s transverse to the heat bath and f t longitudinal to the 
heat bath. The ratio of both determines the pion velocity 



ft' 



(11) 



At temperature T = we have f„ = f s = ft and conse- 
quently u = 1, whereas at nonzero temperature this is no 
longer the case. We stress that the pion decay constants 
constitute static quantities and are in principle calculable 
in a thermodynamic equilibrium approach as pointed out 
in [30l ]. The relation between the pion decay constant, 
the full quark propagator and the pion Bethe-Salpeter 
amplitude in the vacuum is outlined in Ref. [32j ■ To ob- 
tain the decay constants at finite temperatures we may 
proceed along the same lines considering the in-medium 
pion propagation 



D T , 



<Jl + v?(jP- + ml) 



(12) 



This leads to the following expression for the decay con- 
stants 

P u ft = 3tr D Tj2 J -^T„(q,P)S(q + P)^S(q) 

(13) 

valid on the pion mass shell. We used the abbreviation 
Pa — ( u Pj w p)- I n the chiral limit m\ = and the expres- 
sion evaluated in the limit P^ — > determines the decay 
constants: The longitudinal decay component f t is ob- 
tained from the time component of Eq. (|13[) , whereas the 
transverse component / s = uft can be extracted from the 
spatial components. In order to be able to compute them 



4 



and to specify the pion contribution T£ we need the pion 
Bethe-Salpeter amplitude IV It satisfies a homogenous 
Bethe-Salpeter equation which in principle needs to be 
solved at finite temperature. This is a tremendous task 
beyond the scope of the present work. However from the 
axial vector Ward-Takahashi identity in the chiral limit 
it follows that the Bethe-Salpeter amplitude for P 2 = 
may be written as 



r w (9,o) = 75 



B(q) 
ft 



(14) 



This constitutes a Goldberger-Treiman like relation for 
the quark-pion coupling and is known to represent the 
leading term also for P 2 ^ at zero temperature. In the 
following we approximate the Bcthc-Salpctcr amplitude 
for general momenta by 



r w (g,P) «75 



B(g) 
ft 



(15) 



For our purposes this approximation seems particularly 
reasonable since it incorporates the long-range interac- 
tion in Eq. which is important for the critical be- 
havior. Using (fT5|) in (fT3|) yields the Pagel-Stokar ap- 
proximation at finite temperature for the pion decay con- 
stants. 

Following the above considerations and knowledge of 
the gluon and Yang-Mills quark-gluon vertex dressing 
functions provided we obtain a closed system of equa- 
tions for the quark propagator. This truncation also in- 
cludes effects of pion backcoupling. Up to now we did not 
consider the contribution of the scalar flavor singlet me- 
son T°. We know that its correlator is degenerate with 
the 7r-correlator for momenta large compared to its mass. 
Near the transition temperature and close to the chiral 
limit m — > the mass m a tends to zero and therefore 
sufficiently close to T c we expect the above relations to 
hold also for the sigma meson, with 75 in (|15p replaced 
by one. This will be exploited in the following. 



III. ANALYTICAL SCALING ANALYSIS 

Having outlined general properties of pions in the sym- 
metry broken phase we will now perform a scaling analy- 
sis of the quark DSE given in Eq. (TTJ) for T < T c close to 
the critical temperature. We concentrate on the scalar 
dressing function B(p 7 u} p ) and work with very small but 
nonzero bare quark mass m. The reduced temperature 
is given by t=(T c - T)/T c . 

The scaling behavior of the pion velocity and decay 
constants close to the critical temperature can be ob- 
tained from a matching of an effective theory with QCD 
at the scale m a . As shown by Son and Stephanov in 
Ref. this matching reveals 



using the known scaling relations for the inverse correla- 
tion length and the order parameter 



t v 



v/2(l- v ) 



(17) 



In d = 3 dimensions the values for the exponents of the 
0(4)-universality class are v w 0.73 and 77 « 0.03, see 
eg- 10. 

For our scaling analysis we assume the scalar function 
to fulfill a scaling law, B ~ t x , with some unknown expo- 
nent x. On the other hand there is no indication for the 
vector dressing functions A and C to scale with reduced 
temperature. The main idea is to utilize self-consistency 
of the Dyson-Schwinger equation. This requires left- and 
right-hand side of ([T]) to scale with t in the same way. For 
simplicity we ignore the anomalous dimension rj which 
would be only a small quantitative correction anyway. 
Focusing on B it remains to analyze the scalar projec- 
tion of the self-energies given in Eqs. ©, © and (gj: 



B{t) = EymW 



£f(i) + £f(i), 



(18) 



with SyM = 

Tr(E YM )/4 respectively = Tr(S,, CT )/4. 
The important quantities of the analysis are the scalar 
dressing function B and the transverse pion decay con- 
stant f s . In SyM an d the scalar dressing B occurs in 
the non-perturbativc quark propagator and in the Bethe- 
Salpeter amplitude. The pion decay constant f s only oc- 
curs in E^ j0 . in the pion velocity u = f s / ft- We stress 
that the scalar dressing function in the denominator of 
the quark propagator can safely be ignored for our anal- 
ysis, since B is very small close to T c and therefore dom- 
inated by the large, non-vanishing fermionic Matsubara 
frequencies in the denominator. 

For the first contribution we find Sy M (t) ~ B(t), since 
only B(q,(jj q ) in the integral kernel of 5]y M scales with 
reduced temperature 2 . Therefore self-consistency is triv- 
ially fulfilled and no constraint for x can be derived from 
this. On the other hand, the projection of the self- 
energies yields 



B{t) ~ (i) 



B 3 (t) B 5 (t) 



mt) fm 



(19) 



where c is a dimensionful constant with respect to t. The 
numerator is obtained from the fact that the trace over 
Dirac matrices is only non-vanishing for an even number 
of matrices. The denominator stems from the zero fre- 
quency meson propagator with lo 2 = 0. For small B, i.e. 



v/2 



(16) 



2 Here we already use that the Yang-Mills vertex dressing and 
the gluon dressing function are independent of t. The latter 
assumption can be justified by a scaling analysis of the quark- 
loop contribution to the gluon propagator along the lines fol- 
lowed here. One consistently obtains scaling with t proportional 
to B 2 /u 2 ~ const, cf. Eq. lO- For the Yang-Mills part of 
the vertex one can show that close to T c diagrammatic contri- 
butions containing only dressing functions A, C dominate and 
consequently this part of the vertex also does not scale. 
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close to T c , the cubic term is leading. With the scaling 
of the transverse pion decay constant given in Eq. (lit 
we find 



t x 



t v 1 



(20) 



Thus our analysis suggests that in case of a second order 
phase transition of Nf = 2 QCD the quark scalar dress- 
ing function should be found to behave like W 2 with 
v w 0.73. Our result is supported by the fact that this 
solution is not only consistent with the quark propaga- 
tor Dyson-Schwinger equation but also with the known 
relations for the quark condensate (|17[) and for the trans- 
verse pion decay constant (fT5]) . To see this we recall that 
the quark condensate can be obtained from the full quark 
propagator by 



including meson back-coupling effects. This constitutes a 
two loop expression whose numerical computation is still 
non-standard and beyond the scope of this work. There- 
fore and in order to avoid technical problems arising in 
such a computation we simplify the expression further. 
This is done similar to Refs. [20|. We first give the de- 
tails of this approximation and then check the analysis of 
the previous section by setting a power law for the decay 
constant f s . 



FIG. 3: The approximated Schwinger-Dyson equation for the 
quark propagator. 



' N - T J2 j (^3^(p^ p ) 



(21) 



with a B (p,u p ) = B/(oj%C 2 (p,u; p ) + p 2 A 2 (p,uj p ) 
B 2 (p,uj p )). Therefore 



(22) 



consistent with (|17[) in the approximation of vanishing 
anomalous dimension. For the pion decay constant we 
conclude by inserting B ~ t v l 2 in Eq. (|T3|) and evaluation 
for P = (p, 0) 0: 



B ~ t v l 2 => f. 2 ~ t v 



(23) 



Hence, we obtain a closed system in terms of scaling with 
reduced temperature t. 

The above analysis relies on self-consistency of the 
quark-DSE. Given a non-trivial scaling behavior of / s , 
we find scaling of the scalar part of the quark propaga- 
tor, which is induced via the pion dispersion relation. We 
therefore identified the mechanism for non-trivial scaling 
of the quark. However, the emergence of non-mean-field 
critical exponents in f s is not explained. This issue needs 
to be further clarified in future studies. In the following 
we verify the analytical findings of this section also in a 
numerical treatment of the quark-DSE. 



IV. NUMERICAL SOLUTIONS OF THE DSE 

In the previous section we presented a self-consistent 
approximation for the quark Dyson-Schwinger equation 



Truncation scheme 



In order to further simplify the diagram in the lower 
line of Fig. [5] note that one of the loops involves two bare 
quark-gluon vertices, a dressed gluon propagator and a 
full Bcthc-Salpeter vertex. If the gluon propagator and 
the two quark-gluon vertices would be replaced by a full 
Bethe-Salpeter kernel this loop would exactly match a 
Bethe-Salpeter equation. From this we deduce that the 
leading part of this loop is given by a term proportional 
to Z<z^t % in the pseudoscalar meson sector and Z^t 1 for 
a scalar meson. In the vacuum it has been shown in 
Ref. [28| that good results for meson phenomenology can 
already be obtained if the proportionality factor is set to 
one. Here such a simple truncation is not possible, since 
it would destroy the scaling properties close to T c ana- 
lyzed in the previous section. Instead, we need to replace 
the loops by Z2J5T 1 and Z^j 1 , each modified with their 
respective quark-meson couplings r„. j<7 . The resulting ap- 
proximated quark Dyson-Schwinger equation, displayed 
in Fig. [31 has the same scaling properties as its two-loop 
counterpart in Fig. [2] Furthermore it is sufficient for the 
scaling analysis to consider only the zeroth Matsubara 
frequency of the mesons and choose = m a = 0. The 
corresponding equation is then given by 
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/<Pq — 

d 3 q 



-T 



(2tt) 3 



2Ti(p,uj p )S(q,w p )T v (q,u p )D„{k,0)+Ti{p,uj p )S(^^ . (24) 



As an aside, note that this truncation also fulfills 
the Mermin- Wagner theorem: Using Eq. (| 14[) for the 
Bcthc-Salpcter amplitude we find that in case of sponta- 
neous symmetry breaking the Goldstonc boson propaga- 
tor D n (k,0) would lead to an infrared singular integral 
in Eq. (|24p for d < 2 space dimensions. The singularity 
would shift B(p,oj p ) — s- — oo. Therefore in d < 2 space 
dimensions only the symmetric phase with B{p,u> p ) = 
can be realized in agreement with Mermin- Wagner. 

To implement Eq. (|24l) it remains to specify the input 
for the Yang-Mills part of the quark-gluon interaction 
and the gluon propagator. Since we are looking at critical 
scaling related to universality we expect that the details 
of both are not crucial for our studies. We have checked 
this by using two different approaches: On the one hand 
we used a simple rainbow-ladder like phenomenological 
model, which is detailed in appendix |XJ On the other 
hand we followed Ref. [22[ and used quenched lattice data 
for the temperature dependent gluon propagator together 
with an explicit back-coupling of the quarks onto the 
gluon propagator as input, see appendix [A] for more de- 
tails. Indeed we have found that the results for the crit- 
ical scaling behavior close to T c are similar in both ap- 
proaches, as expected. In the next section we therefore 
concentrate on one of the two approximation schemes and 
present results for the more advanced truncation with 
temperature dependent gluon only. 



Results: scaling of quark dressing and 
condensate 



In order to study scaling close to the critical temper- 
ature we first have to determine T c . This can be done 
either from the behavior of the scalar quark dressing func- 
tion B or, equivalently, from the chiral condensate. Our 
results for the former are shown in Fig. [4] We observe 
a reduction for B(p = 0,u> p = irT,T) with temperature 
which results in a phase transition at T c = 162.5 ± 2.5 
MeV, where the error reflect the coarseness of our tem- 
perature grid. For small temperatures the details of this 
reduction do depend on the details of the interaction, 
since the Yang-Mills part of the quark-gluon interaction 
is leading there. For temperatures close to the critical 
one, effects from meson back-coupling become more and 
more important due to the long range nature of the me- 
son dispersion relations. In Fig. U we have been using 
all Matsubara frequencies in the pion propagators. Ne- 




0.1 
T [GeV] 



FIG. 4: The variation of the scalar quark dressing function 
B(p = 0,w. 

chiral critical temperature T l 



p — tyT) with temperature not too close to the 



glecting all but the zeroth Matsubara frequency (which 
is sufficient for scaling around T c ) and considering also 
the sigma propagator which becomes important close to 
T c we obtain a similar picture with only slightly changed 
temperature T c = 157 ± 2.5 MeV. We use this value for 
the following scaling analysis. 

As shown analytically in section UTTI we obtain the cor- 
rect 0(4)-scaling in the quark-DSE provided we feed in 
the correct scaling properties of the spatial part of the 
pion decay constant f s . In order to numerically test this 
scenario we choose the following scaling ansatz for f s 



fs(T) 



/.(To) 



(T c - To)-/ 2 



(T c - T) 



v/2 



(25) 



where T = 150 MeV and f s (T ) is the result for the 
decay constant obtained from Eq. (| 13[) . The critical tem- 
perature is chosen to be T c = 162 MeV. With this con- 
struction our numerical results smoothly connect with 
the ones for f s (T) at T < 150 MeV. The exponent is 
set to v = 0.73 corresponding to the value of the 0(4) 
universality class. 

In figure [5] the logarithm of the zeroth mode of the 
scalar dressing function B (red circles) at zero momen- 
tum is shown as a function of ln(tT c ). In addition, the 
proportionality to t v l 2 is shown by the solid line. For 
temperatures close enough to the critical one, B{p = 
0, Lo p — 7rT) behaves exactly as it should according to our 
analytical scaling analysis in section IIII1 We also show 
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FIG. 5: Scaling of the scalar quark dressing function B(p = 
0, cj p = 7vT) as well as the condensate with reduced tempera- 
ture t = (T c - T)/T c . 



results for the chiral condensate (blue squares) calculated 
from Eq. (j!?Tj) . They likewise follow such a scaling behav- 
ior as seen by the straight line ~ W 2 . The size of the 
scaling region can be inferred from the deviations from 
scaling at large t which set in at roughly T c — T « lOMeV. 

Finally we investigated, whether the scaling of f s can 
be obtained from its defining equation (|13[) . To this end 
we abandon the scaling ansatz but calculate f a self- 
consistently from Eq. (fT3|) . It turns out, that in this case 
we do observe scaling around T c , but the corresponding 
anomalous dimensions are the mean field ones. We at- 
tribute this result to the approximations done in the pion 
sector, which reduce Eq. ([Tl?]) to the Pagels-Stokar form 
which may not be sufficiently accurate close to T c . 



V. CONCLUSIONS AND OUTLOOK 

In this paper we investigated the quark propagator 
close to the second order chiral phase transition for 
Nf = 2 massless quarks in the symmetry broken phase. 
We used the continuum formulation of the theory utiliz- 
ing the Dyson-Schwinger functional framework. In con- 
trast to other methods we have the ability to work in the 
chiral and thermodynamic limit within this approach. In 
order to take into account the relevant degrees of freedom 
we outlined a truncation scheme for the quark propaga- 
tor including meson backcoupling effects. In particular 
we studied the effect of the Goldstone bosons, the pions, 
and the radial excitation, the sigma meson. 

Identifying the quantities of our truncation obeying 
known scaling laws we analytically derived the scaling 
law of the scalar dressing function B with reduced tem- 
perature. Our analysis reveals that the behavior of the 
pion velocity u respectively the pion decay constant f s 
is crucial for the scaling of the scalar dressing function 



close to T c . Furthermore within the employed approx- 
imation (neglecting rf) the result is consistent with the 
known scaling relations for the chiral condensate and the 
decay constant. 

In addition we also considered the quark Dyson- 
Schwinger equation numerically. Here we could verify our 
analytical findings also numerically: using an appropriate 
scaling ansatz for the pion decay constant we showed that 
the scalar quark dressing function and the chiral conden- 
sate scale with 0(4)-critical exponents in a temperature 
region of about 10 MeV below the chiral critical tem- 
perature. Using two different truncation schemes for the 
Yang-Mills part of the interaction we explicitly verified, 
that chiral critical scaling is independent of the details 
of the gluon propagator and the Yang-Mills part of the 
quark-gluon vertex, as expected from universality. 

Finally we performed a first attempt to self- 
consistcntly close the system by determining the pion de- 
cay constant from the Pagel-Stokar formula. This, how- 
ever, turned out not to be sufficient, since the resulting 
critical exponents returned to their mean field values. We 
interpret this a a failure of the Pagel-Stokar approxima- 
tion close to the chiral phase transition. Going beyond 
this approximation will be a task for future work. 
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Appendix A: Two truncations of the Yang-Mills 
part of the quark-DSE 



In the main part of this work we numerically inves- 
tigate a combination of purely gluonic effects and me- 
son back-coupling effects onto the quark propagator. 
Since the Yang-Mills part of the quark-gluon interac- 
tion does not take part in the critical 0(4)-scaling close 
to T c , it is sufficient for our purpose to approximate 
it in the simple fashion of a rainbow-ladder like model 
[35| . For the Yang-Mills part of the quark-gluon vertex 
rJ M = 7^r YM (z) with the squared gluon momentum 

z = k 2 = ul + k 2 = (uj q - Lo p ) 2 + (q — p) 2 we use the 



8 



following model ansatz 



r YM (z) = 



z \ I «1 



z + d 2 



zd 3 



d 2 + z dl + {z-d 2 ) 



+ 



^QCD + z 



4tt 



A 2 

QCD 



/W \\og (z/A 



QCD 



QCD 



-2(5 



(Al) 



with parameters g?i = 1.45GeV 2 , d 2 = O.lGeV 2 and d$ = 
3.95GeV 2 . Furthermore the QCD scale is A QC d = 0.520 
GeV, = 0.2 and k « 0.595. Determined from per- 
turbation theory we find 5 = ~9N C /(4:4:N C - 8Nf) and 
/?o = (HJVc — 2Nf)/3. The gluon propagator dressing 
functions are approximated by a fit function to the zero 
temperature results 

Z(z) = Z T {z) = Z L (z) 



Z(z) = 



4tt 



2 + ^QCD 



z/A 2 



QCD 



Po z + A 2 C 



QCD 



ln(z/A 2 



QCD) 



5 Mc 



- 1 



(A2) 



where /3q = (11AT C — 2Nf)/3. The different temperature 
dependence of the longitudinal and transverse dressing 
functions is neglected. 



As an alternative we also employ the much more com- 
plex interaction developed in [22j. There, quenched lat- 
tice data for the temperature dependent gluon propaga- 
tor have been combined with the quark loop polariza- 
tion diagram to obtain a numerical expression for the 
unquenched gluon propagator with Nf = 2 quarks back- 
coupled. Unfortunately, there is no reliable information 
on the details of the quark-gluon vertex at finite tem- 
perature, so also in this approach one has to resort to a 
model for the vertex. Here we choose 



FIG. 6: The DSE for the gluon propagator. Filled circles 
denote dressed propagators and vertices. The shaded circle 
denotes the quenched gluon. 



T u (q, k,p) = Z 3 j u 



di 



A 2 



/3 a(//)ln[g 2 /A 2 + l] 



4tt 



2ri\ 



(A3) 



r 



where q = (q,uj q ) denotes the gluon momentum and p = 
(p,uip), k = (k,Uk) the quark and antiquark momenta, 
respectively. Furthermore 26 = -18iV c /(44JV c - 8N f ) is 
the anomalous dimension of the vertex and the parame- 
ters arc di = 15GcV 2 ,d 2 = 0.5 GeV 2 , A = 1.4 GeV and 
we renormalize at a(fi) = 0.3. 

The numerical results in section ITVB I of this work have 
been obtained using the second, more complicated trun- 
cation of the Yang-Mills part of the quark-gluon interac- 
tion. However, we wish to emphasize once more that the 
results using the simpler rainbow-ladder form, (|Aip and 
(|A2[) . are similar. Critical scaling of the chiral order pa- 
rameters is universal and thus independent of the details 



of the Yang-Mills part of the quark-gluon interaction. 
In particular this is true for the deep infrared behav- 
ior: whereas the rainbow ladder approximation (|A1[) and 
(|A2[) is constructed to adhere to the Yang-Mills 'scaling' 
scenario (in the terms of Rcf. [36]), the more compli- 
cated construction (|A3[) together with the t emp erature 
dependent lattice gluon propagator of Ref. [22( behave 
according to the 'decoupling' scenario (again in terms of 
Ref. [36[). Although the deep infrared behavior of these 
two scenarios is very different, it does not affect the chiral 
critical scaling properties of the systems at T c , in agree- 
ment with previous findings (2lj . 
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